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Abstract 

We provide a detailed analysis of the supersymmetric tri-lepton signals for spar- 
ticle searches at the Tevatron in the minimal supersymmetric standard model with 
general CP-violating phases but without flavor mixing among sfermions of different 
generation. The stringent experimental constraints on the CP-violating phases from 
the electron and neutron electric dipole moments are included in the analysis for two 
exemplary scenarios of the SUSY parameters; one with decoupled first two genera- 
tion sfermions and the other with non-decoupled sfermions. In both scenarios, the 
production cross section and the branching fractions of the leptonic chargino and 
neutralino decays are sensitive to CP-violating phases. The production-decay spin 
correlations lead to several non-trivial CP-even observables such as the lepton in- 
variant mass distribution and the lepton angular distribution, and several interesting 
T-odd (CP-odd) momentum triple products. The possibility of measuring the CP- 
violating phases directly through those T-odd observables is investigated in detail. 
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1 Introduction 



Supersymmetry is the currently best motivated extension of the Standard Model (SM) of 
particle physics which allows to stabilize the gauge hierarchy without getting into conflict 
with electroweak precision data. Among all possible supersymmetric theories, the Minimal 
Supersymmetric Standard Model (MSSM) occupies a special position. It is not only the 
simplest, i.e. most economical, potentially realistic supersymmetric field theory, but it also 
has just the right particle content to allow for the unification of all gauge interactions 

The R-parity preserving MSSM as a softly-broken SUSY model contains in general 
more than one-hundred independent physical parameters including about forty-five CP- 
violating physical phases in the Lagrangian. The CP-violating phases, if they are large, 
can give a significant impact on not only various CP-odd observables but also the CP-even 
quantities such as sparticle masses [Q], production cross sections, branching fractions [Q, 
LSP relic density 0, Higgs boson masses and couplings 0, CP violation in the B and K 
systems [0, and so on. 

The most stringent (indirect) constraints ^ on the MSSM CP-violating phases come 
from the precise measurements of the electron and neutron electric dipole moments (EDM), 
but the indirect EDM constraints depend strongly on the assumptions taken in the analy- 
sis without any a priori justifications. Relaxing the rather stringent assumptions, several 
recent works [p!0| , pj| have shown that the constraints could be evaded without suppressing 
the CP-violating phases of the theory. One option |11| is to make the first two genera- 



tions of scalar fermions very heavy so that one-loop EDM constraints are automatically 
evaded while keeping the third-generation sfermions relatively light to preserve naturalness. 



This case can be naturally explained by the so-called effective SUSY models |]T2[ where de- 
couplings of the first and second generation sfermions are invoked to solve the SUSY FCNC 
and CP problems without spoiling naturalness. Another possibility is that various SUSY 
parameters are arranged ||10| to lead to partial cancellations among their contributions to 
the electron and neutron EDMs without taking very large sfermion masses. Consequently, 
it is not yet clear at all whether the CP-violating phases of the MSSM are small or not. 
Once supersymmetric particles are discovered at colliders, it will be therefore of great im- 
portance to directly measure the CP-violating phases as well as the other real parameters 
of the supersymmetric Lagrangian. 



In this paper we illustrate how the presence of the CP-violating phases affects observ- 
ables which can be measured in the classical reaction pp — >• xf X2 with subsequent decays 
Xi — >■ Xi^^^e X2 ~^ Xi^'~^^'~ {^1 ^ = e, n) when the lightest neutralino Xi is the 
lightest supersymmetric particle (LSP). Such tri-lepton signatures have been investigated 
in several CDF and DO analyses at the Tevatron based on the models with only real SUSY 
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parameters. That is to say, most of the works [O have been done under the assumption 
that all the couplings are related at the grand unification or Planck scale and they are real. 
In the light of the possibility of large CP-violating phases, we investigate systematically 
the impact of the CP-violating phases on the SUSY tri-lepton signatures at the Tevatron, 
including the constraints on the phases imposed by the electron and neutron EDMs and 
taking into account the full spin/angular correlations between the production and the lep- 
tonic decays of the associated chargino and neutralino pair. 



The paper is organized as follows. In Sect. we describe the chargino, neutralino and 
flavor-preserving sfermion mixing phenomena on the most general footing and identify all 
the relevant physical CP-violating phases. Section |^ is devoted to the detailed discussion 
of the constraints by the electron and neutron EDM measurements on the CP-violating 
phases. In Sect. ^ we present in detail the formalism to describe the associated produc- 
tion of the chargino xf cind the neutralino X2 PP collisions at the Tevatron; production 
helicity amplitudes, chargino and neutralino mass spectra, and chargino and neutralino 
polarization vectors. Section || is devoted to the detailed description of the decay modes 
of the chargino xf the neutralino X2 their branching fractions. In Sect. ^ after 
explaining the method of obtaining the fully spin-correlated distributions of the associated 
production and decays of the chargino and neutralino, we investigate in detail the impact 
of the CP-violating phases on various physical observables such as the total rate of the tri- 
lepton signatures, the dilepton invariant mass distributions, the lepton angle distribution 
as well as the CP-odd triple products of three proton and/or lepton momenta. Finally, we 
summarize our findings and conclude in Sect. |^. 



2 Supersymmetric Flavor Conserving Mixing 

The existence of the non-trivial MSSM CP-violating phases is due to SUSY breakdown, 
so that the CP-violating phases appear in the soft-breaking parameters and the mixing 
among sparticles due to the electroweak gauge symmetry breaking. As a whole, the MSSM 
has three well-known sources of CP violation. The first is related to the two Higgs-boson 
doublets present in the model as both the /i parameter in the superpotential and the soft 
breaking parameter B can be complex. We denote the phases of /i and B by $^ and 
respectively. Secondly, there are three more phases {^1,^2,^3} related to the com- 
plex U(l)y, SU(2)i and SU(3)c gauginos masses. Finally, most of the other CP-violating 
phases originate from the flavor sector of the MSSM Lagrangian, either in the scalar soft 
mass matrices or in the trilinear matrices. 

The sfermion mass matrices are Hermitian so that only off-diagonal terms can be com- 
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plex, but the trilinear matrices are in general 3x3 matrices allowing for the complex 
diagonal entries. The effects of the phases associated with the off-diagonal terms on exper- 
imental observables are strongly suppressed by the same mechanism required to suppress 
the flavor changing neutral current effects. Therefore, we neglect these flavor-changing 
CP-violating phases in the present work assuming that all the scalar soft mass matrices 
and trilinear parameters are flavor diagonal and the complex trilinear terms Af with its 
phase are proportional to the corresponding fermion Yukawa couplings. Consequently, 
the gaugino masses Mi, M2, M3 and the higgsino mass parameter n as well as the trilinear 
parameters Af can be complex in the CP-noninvariant theories. However, reparameterising 
the flelds one can take M2 to be real and positive without any loss of generality and all other 
parameter choices are related to the speciflc choice by an appropriate R transformation. 

Neglecting flavor mixing among sfermions, the sfermion mass matrix squared is given 

by 



M 



m 



fLR ^f^+^f + Dfu 



where the mixing term rn^j^j^, and Df^ and D^^ are: 

rnjj^j^ — —m f (^A*f + Ijl tan (3/ cot (3^ for f — d,e/u, 

Df^=m\cos2p{Ti^-Qfsl) , D = m| cos 2/3 . (2) 

The flrst term of each diagonal element is the soft scalar mass term evaluated at the 
weak scale, and the second is the mass squared of the corresponding sfermion mass, and 
the last one is the so-called D term of the MSSM superpotential. The trilinear term Af 
causing the left-right mixing is due to the soft-breaking Yukawa-type interaction, and ^ 
is the complex supersymmetric higgsino mass parameter and tan/? is the ratio V2/V1 of the 
vacuum expectation values of the two neutral Higgs fields which break the electroweak gauge 
symmetry. The sfermion mass eigenvalues and eigenstates can be obtained by diagonalizing 
the above mass matrix with a unitary matrix Uf such that UjjA^jUf = diag(mj^, m^^). 
We parameterize Uj as 

^ f cos Of - sin Ofe-^t'f \ 
f \ sine fc^'t'i cosOf I' ^' 

where 0/ = arg[— m/ {Af + //*tan/3)], < ^/ < tt and — 7r/2 < (t)f < 7r/2 without any loss 
of generality. 
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In the MSSM, the spin-1/2 supersymmetric partners of the boson and the charged 



Higgs boson, W"^ and H"^, respectively, mix to form the chargino mass eigenstates xt 
The chargino mass matrix dictating this mixing is given in the H~} basis by 



M 



c 



M2 \/2mwCp 



\/2mwSi3 |//|e' 



(4) 



which is built up by the fundamental SUSY parameters; the SU(2)/, gaugino mass M2, the 
modulus and phase of /x, s^j = sin (5 and cp = cos p. Since the chargino mass matrix Mc is 
not symmetric, two different unitary matrices acting on the left- and right-chiral {W, H) 
states are needed to diagonalize the mass matrix Aic'- 

(5) 

L,R \ / L,R 

SO that UrM-cUj^ — diag (m-±, m~±) with the ordering of m~± < m~± as a convention. 





X2 



X2 



The supersymmetric spin-1/2 partners of the neutral U(l)y and SU(2)x, gauge bosons, 
B and W^, respectively, mix with the supersymmetric fermionic partners of the neutral 
Higgs bosons, and to form mass eigenstates. The four physical mass eigenstates 
Xi (i = 1 to 4), called neutrahnos, are obtained by diagonahzing the 4x4 neutrahno mass 
matrix 



M 



N 



I 



\ 



|Mi|e^*i 


-mzcpsw 
mzspsw 





M2 

mzCfsSw 
-mzspcw 



-mzcpsw 
mzcpcw 




mzSjjSw 

-mz-SfsCw 





(6) 



where sw = sin 9w and cw = cos6'^y. Since the neutralino mass matrix A^at is a complex 
and symmetric so that it can be diagonalized by a single unitary matrix N as N*Ai AriV" = 
diag (m^j-o, m^o, rrij^o, m^o) with the ordering of m^o < m^o < m^o < m^o. On the other 
hand, the gluinos, the spin-1/2 partners of gluons, do not mix among themselves and in- 
volve only a single complex phase $3. 



To recapitulate, we have in total twelve non-trivial CP-violating phases; two phases 
{$1, $3} from the gaugino sector, one phase $^ from the higgsino sector, and nine phases 
for the nine sfermion flavors. 
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3 Electron and Neutron Electric Dipole Moments 



The electric dipole interaction of a spin-1/2 charged particle / with an electromagnetic 
field is described in a model-independent way by the 5-dimensional effective Lagrangian 

Cedm = -'-dfla'^'j^fF,,. (7) 

In theories with CP-violating interactions, the electric dipole moment df receives con- 
tributions from loop diagrams. In the case of the electron EDM the chargino-sneutrino 
and neutralino-slepton loops contribute whereas in the quark case the chargino-squark, 
neutralino-squark and gluino-squark loops are involved (See Fig. |l|). Recently, one-loop 
chromoelectric dipole moment (CEDM) contributions have been studied in Ref. and 
found to be comparable to the EDM ones. In addition, there is a significant contribution 
to the neutron EDM through the so-called Weinberg's three-gluon dimension-six operator 
|1T0| , . Therefore, for the CEDMs of quarks we include the chargino-squark, neutralino- 
squark and gluino-squark loops, where the gluonic dimension-six operator gets contribu- 
tions from the loop containing top (s)quark and gluinos. The quark CEDM is then defined 
as the coefficient dg in the 5-dimensional effective Lagrangian 

/:cEDM = 4 ga'-'l^T'^q , (8) 

where the indices a (= 1 to 8) denote the gluons and T°- are the SU(3)c generators. And, 
the gluonic dimension-six operator is given by 

Ag = ~ g dcjahc ^^"^^ G'^^pG^SG'i^ . (9) 

The various SUSY contributions to the fermion EDM are illustrated in Fig. ^ the diagram 
(a) denotes the chargino, neutralino, and gluino loops where the dashed line of each loop 
is for a scalar fermion mass eigenstate corresponding to the fermion, and the diagram (b) 
is for the effective 3-gluon operator. Both of them contribute to the coefficients df and dg 
depending on whether a photon or a gluon is radiated off the loop. These contributions can 
be calculated at the electroweak scale since a typical SUSY scale in most SUSY models is 
of the same order of magnitude as the electroweak scale. 



In the following we enumerate all the SUSY contributions to the first -generation fermion 
EDMs. The chargino-loop contributions are determined by the Lagrangian 



Sw 



1=1 j=l 



(10) 
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describing the X^^f^f interactions for the first generation sfermions without flavor mixing 
Here Pl,r = (1 T 75)/2 and {Yg, Y^, Yy} are the Yukawa couphngs: 

Trip rrid m. 



Y.^ 



\/2mw S/3 



(11) 



One can readily find out that the interactions of the charginos with a fermion and a 
sfermion contribute to the EDM of the fermion / as 

,2 



ea 



Atts 



EE 



2 

A, 



3 Cj! 



2 

W i=i j=i "^J. 



m 



QfB\ 



m 



fi 



+ {Qf-Qf)A 




where the dimensionless functions A{r) and B{r) are defined as 



2(1 -r)2 



3-r + 2 



Inr 
1 — r 



B{r) = 



2(1 -r)^ 



l + r + 2 



r Inr 
1 — r 



(12) 



(13) 



and 



Af" ^Y, U 



Ulj2UlijiSii 

Yu Un.^ui [Ul^.ui* - Ya Ul^.ut 

TT* TTU* V TT* TJU"* 



{Qf = Qe, Qf = 0, TTlf. = mp^), 

{Qf = Qu, Qf = Qd, = rriui), 
{Qf = Qd, Qf = <5«, ruf. = mj.). 



(14) 



On the other hand, the chargino contributions to the CEDM of the quark q are given by 



9sOi 



2 2 777, 



w 



/ 2 



-ij=i 



2 



(15) 



Secondly, the neutralino-loop contributions to the EDM of a fermion / are determined 
by the Lagrangian describing the most general x°-/-/ interactions, 

2 4 



^x°// = "7^ EE/ + B!-Pn] xYh + H.C., 

\/ ISw i=i j=i 

where the couplings B-j' and B-j^ are given by 

4^ = V2YfN;^Uf,^ - 2QfN;, tan ^^C/;^, , 

TiNj2 + {Qf- Ti) iBaiOwNji] Uf^, + v^YfNjhUf^, , 



(16) 



4" = 2 



(17) 



with h — 3 ior f — d,e and h — A ior f — u, respectively. The neutralino contributions to 
the EDM of a fermion / is then given by 



d 



ea 



-Q/EE -*»nr 




(18) 
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and the neutralino contributions to the CEDM of the quark q can be obtained by simply 
replacing e by gg and Q/ by unity in eq. (|^). 

Thirdly, the gluino-loop contributions to the fermion EDMs are determined by the 
Lagrangian describing the most general g-q-q interactions 

C-,,, = V2g,{T'^)i^ (e-^"-^/' iPLTqTR + e'"''' iPRr~qTi) + H.c. , (19) 



where the indices / and m (= 1 to 3) are for the color of quarks and squarks, the index a 
(= 1 to 8) is for the color of a gluino and the 3x3 matrix is a SU(3)c generator. The 
gluino contributions to the coefficients dq and dq are given by 



Opn, 2 m"^ 



4 = ^ E ^^ni(e-^*3^,2.^f .) C{x) , (20) 



i=l ""qi 

with the dimensionless function C{x) defined as 



Bir) 

C(r) = -3A(r) + ^ . (21) 



Finally, the leading nontrivial MSSM contribution to the CP-odd three-gluon term dc 
is given by a two-loop diagram involving the top, top squarks and gluinos pO| : 

3a^a'^ iTT'tinT-'r — ni'^ ) r r 

do = '^' '^m{e-^''^Ul,Ul;)H{x,,X2,x,), (22) 

lovr"^ m~ 

where Xi = mf^/m'^, X2 = rrif^/m'-, x^ = 171^2 /m'~, and the two-loop function H is given by 
the three-fold integral 

1 /"i A^iA^2 

H{zi,Z2,Z3) = - dx du dyx{l-x)u——, (23) 

2 JO Jo Jo 

where the explicit forms of the Ni, N2 and D are given by 

A^i = ^(1 — x) + Z3x{l — x){l — u) — 2ux[ziy + ^2(1 — y)] , 

N2 = {l~xf{l-uf + u^-lx\l-uf, 

9 

D = u{l — x) + Z3x{l — x){l — u) + ux[ziy + Z2{1 ~ y)] ■ 
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Having defined the contributions from tlie individual Feynman diagrams, tlie total EDM 
of electron is given simply by the sum of the chargino and neutralino contributions: 

d, = df + df . (24) 

In principle, the Kobayashi-Maskawa (KM) CP phase in the SM can contribute to the 
electron EDM, but it turns out to be effective only at the three-loop level so that the 
contribution is too small to be considered. 



On the other hand, the EDM of the neutron, a composite particles of quarks and 
gluons, can be obtained from the EDMs of the constituent quarks and gluons. The quark 
EDM evaluated at the electro-weak scale must be evolved down to the hadronic scale via 
renormalization group equation (RGE). Generally, the prescription based on the effective 
chiral quark theory given in Ref. [|15|, are used to calculate the neutron EDM from the 
quark EDMs. According to the model, the neutron EDM is given by 

4 1 

dn = -dd- - du. (25) 

In order to evaluate the numerical value of each individual quark EDM from the effects 
of the EDM operators given in eqs. ([TsD , ([Tsl), (|20|) and (^2]) we use the so-called "naive 



dimensional analysis" as proposed in Ref. [1^. The explicit form of the quark EDMs are 
then given by the contributions of all quark and gluon operators (to leading order in as) 
with the proper dimensional re-scalings as 



d, = {df + d^ + dl) + r/^ - {df + rf^ + rf^) + , (26) 

where 77^, and t]^ are the QCD correction factors due to the renormalization group 
equations, whereas A^^ is the scale of chiral symmetry breaking in QCD; we use t]^ = 1.53 
11, r/^ ~ ^ 3 4 and AsB ^ 1-19 GeV pi. 



Although flavor mixing is neglected, there are still sixteen real parameters and seven 
CP-violating phases: {tan/3, |Mi|, M2, IM3I, mj^, mj^, and {$1, $3, 

for / = e,u,d,t where the SU(2) relation m^^ = m^^ is taken into account. So, any 
reasonable quantitative understanding of the general features of the SUSY contributions 
to the electron and neutron EDMs requires us to make a few appropriate assumptions on 
the SUSY parameters without spoiling their qualitative aspects; (i) we take a universal 
soft-breaking selectron mass rrii and squark mass rriq for the first-generation left- and 
right-handed selectrons and squarks, respectively, and for the left- and right-handed 
top squarks; (ii) the gaugino mass unification condition is assumed only for the modulus 



9 



of the U(l)y and S'U{2)l gaugino mass parameters, i.e. |Mi| = ^taii^ 9w M2 ~ 0.5 M2 and 
the absolute value of the gluino mass, IM3I, is taken to be 500 GeV; (iii) we take \ Af \ to be 
1 TeV for any sfermion / such that the EDM constraints are satisfied. We should be very 
careful in choosing the value of tan f3 which has very significant effect to many SUSY pro- 
cesses. According to the recent calculation for the Barr-Zee-type two-loop contributions 
to the fermion EDMs by Chang, Keung and Pilaftsis [0, the bottom squark contributions 



are very much enhanced for large tan /? so that their contributions cannot be neglected. 
Therefore, for large tan f3 we are forced to introduce additional CP-violating phases in our 
analysis related with the sparticles of the third generation. Moreover, a large tan (3 yields 
a large tau slepton left-right mixing and the Higgs-exchange diagrams in the decays of 
the neutralinos into tau or bottom pairs cannot be neglected due to the enhanced Yukawa 
couplings of third generation. On the contrary, the value of tan/5 less than about 2 has 
been already ruled out by negative results in the Higgs search experiments at LEP II |T9[ 
based on the minimal supergravity scenario with real couplings. However, this experimental 
constraint on tan (3 may be loosened in CP noninvariant theories. Putting off the detailed 
investigation related with the tan (3 dependence of the EDMs, the associated production of 
chargino and neutralino, and their branching ratios, we simply take tan /3 = 3 in the present 
analysis, for which the tau-slepton contributions can be treated on the same footing as the 
other slepton contributions. 

With these assumptions on the SUSY parameters given in the previous paragraph, the 
electron and neutron EDMs are determined by just five real parameters and seven remaining 
CP-violating phases: 

{M2, mg. Trig, mi ; $1, $3, , (27) 

with q = u, d. Let us take two extreme scenarios for allowing relatively large CP-violating 
phases without violating the EDM constraints. Firstly, based on the so-called effective 



SUSY model |]T2[ we decouple the first and second generation sfermions by rendering them 
very heavy without violating naturalness by maintaining the third generation sfermions to 
be light and taking M2 and \fi\ relatively small: 

SI : {M2 = 100 GeV, = 200 GeV, mj = 10 TeV, = 500 GeV} , (28) 

In this scenario, the sfermion masses mj (/ = e, u, d) are large enough to suppress the 
sfermion contributions to the electron and neutron EDMs completely. As a result, the 



present EDM bounds ||20| on the electron and neutron EDMs 

|4| < 4.3 X 10"^^e ■ cm, |c/„| < 1.1 x 10"^^ e ■ cm , (29) 

put no constraints at all on the CP-violating phases {$1.3, "I*^}. Secondly, for the sake 
of generality, we consider a scenario with relatively small universal soft-breaking sfermion 
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masses but with a large value of \fi\ as 

S2: {M2 = 100 GeV, = 700 GeV, mg = 200 GeV, m-^^j = 500 GeV}. (30) 

This scenario of small sfermion masses requires some cancellations among the different 
contributions so as to escape the electron and neutron EDM constraints. The degree of 
cancellations depends on the moduli of the SUSY parameters such as \Mi\ and M2. 
The reason for taking a large (= 700 GeV) is because it allows a relatively large region 
for the CP phases $^ and $1 for any values of the phase between to 27r. The fact 
that the allowed region of two phases increases with has been pointed out in several 
previous works |jlO[- In Fig. 0(a), the allowed range of $^ versus at 95% confidence level 



is displayed with the other real parameters of S2 and the other phases sampled randomly 
within their allowed ranges. The overall trend is that for larger |/i| it is much easier to 
satisfy the EDM limits and any \^\ > 650 GeV allows the full range of Except for small 
the neutron EDM constraints are less stringent than that of the electron. Fig. 0(b) 
shows the allowed region at 95% confidence level for the phases and $1 in the scenario 
S2. Note that the electron EDM constraints allow the full range of $^ only for $1 around vr 
while the neutron EDM constraints give no restrictions to $^ and $1 in the scenario S2. It 
is also worthwhile to note that the allowed region of two phases can be enlarged by taking 
large values of M2 while keeping relatively small. From the above discussion, it is clear 
that allowing large CP violating phases against the stringent electron and neutron EDM 
constraints needs some sort of "fine tuning" among the SUSY parameters. 



On the other hand, the spin-correlated production and decays of the associated chargino 
and neutralino pair in pp collisions is mainly dependent on the flavor-independent CP- 
violating phases and $1 but nearly independent of the flavor-dependent phases 
and the phase $3 of the gluino mass. In the following numerical analysis. Fig. |](b) will be 
taken to be the basic platform for all the contour plots of the production and total cross 
sections of the correlated process, the branching fractions of the chargino and neutralino 
decays, and the CP-odd (T-odd) triple momentum products in the scenario S2. Since the 
neutron EDM constraints on the CP-violating phases are weaker than the electron EDM 
ones in both scenarios SI and S2 we include only the latter constraints on the CP- violating 
phases in our numerical analysis. 



4 Production of a Chargino and Neutralino Pair 

The production of xfx2) followed by the subsequent decays xf ~^ Xi^^^i S'-nd X2 ~^ 
Xii'~^i'~, is a main source of three charged leptons (e or fi) and in PP collisions, called 
tri-lepton events. Since the tri-lepton signal suffers by very tiny SM backgrounds, it is, 
therefore, considered to be one of the most promising channels by which low energy SUSY 
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can be discovered in hadron colliders, especially at the Tevatron. The cross section of 
X1X2 determined mainly by the chargino and neutralino masses and the overall effi- 
ciency for detecting the three final-state leptons is very sensitive to the mass splittings 
between {xf, X2} Xi &s well as the corresponding branching fractions. In the following 
subsection, we investigate the dependence of chargino and neutralino mass spectra on the 
CP-violating phases $^ and $1 in the two scenarios SI and 52. Then, we present the 
helicity amplitudes for the parton-level process du allowing us to construct the 

production cross section and to investigate the chargino and neutralino polarizations. 



4.1 Chargino and neutralino masses 

The chargino and neutralino masses are obtained by diagonalising the mass matrices A4c 
and Mn, respectively. Unlike the neutralino masses, the chargino masses are independent 
of the modulus and phase of the U(l)y gaugino mass, \Mi\ and $1. Since, in both scenarios 
SI and S2, the higgsino mass parameter \fi\ is larger than the gaugino mass parameters 
M2 and I Ml I, the lightest mass eigenstates are gaugino-like and the heaviest states are 
higgsino-like. This feature is expected to be more prominent in the scenario S2 with the 
larger value of |/^|. 

Figure ^ shows the mass spectrum of the lightest chargino xf and the neutralinos x? 2 
on the {$^, $1} plane for the scenarios SI and S2. Each shaded area of the lower figures is 
excluded by the electron EDM constraints. Except for the region of $^ = 0, 27r, the masses 
m-± and m^o are very similar in size and independent of $1 in both SI and S2 while m^o 
exhibits a strongly correlated dependence on the CP-violating phases. The masses m-± 

Xi 

and m^o increase as $^ approaches vr, while m^o becomes maximal at certain non-trivial 
values of $^ and $1 in SI. This implies that m^o is strongly affected by a small value 
of while m-± and m^o are essentially determined by the SU(2)l gaugino mass M2. 
The mass m^o becomes smaller as the CP phases $1 and $^ approach the off-diagonal 
line on the plane, implying that the mass is a function of the sum $^ + $1 of two CP 
phases to a very good approximation. Note that the masses m -± and m^o ^ are more sen- 
sitive to the phases in SI because is comparable with M2 and mz in size in this scenario. 



Although the chargino and neutralino masses cannot be separately measured at the 
Tevatron, the mass difference = rriy^o — m^o can be determined with a good precision 
2T[] by measuring the end points of the invariant mass distribution of the same flavor 



but opposite sign dileptons from the neutralino decay. The precision of determining the 
mass difference is sensitive to the event rate and dependent on experimental capabilities. 
Deferring the discussion on these aspects to Sect. ^ we illustrate what information on the 
CP-violating phases the A^, measurement can provide us with. For that purpose, we 
assume the uncertainty for the mass difference to be 400 MeV and exhibit in Fig. ^ the 
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allowed area of the CP-violating phases {$^, $1} for the SUSY parameters of (a) SI and 
(b) iS2. The figure clearly shows that the mass difference is a very sensitive probe of the 
CP-violating phases if the other real parameters are known; the sensitivity of the mass 
difference to the phases are enhanced when the gaugino and higgsino mass parameters are 
comparable to rriz in size. 



4.2 Parton— level production helicity amplitudes 

Although we are mainly interested in the production process du X1X2' discuss in this 

section the associated production of any chargino and neutralino pair du XiX^j = !> 2 
and j = 1 to 4) for the sake of generality. 

The parton-level production process du xix'j is generated by the three mechanisms 
shown in Fig. ^ the s-channel W~ exchange, the t-channel d^ exchange and the w-channel 
ul exchange. Before presenting the explicit form for the production helicity amplitudes, we 
note that the chirality mixing of the first and second generation sfermions is proportional 
to the fermion mass much smaller than the beam energy of the Tevatron. Therefore, we can 
safely ignore the sfermion left-right chirality mixing in calculating the associated chargino- 
neutralino production rate so that the trilinear term Af does not play any role in the 
high energy process. The phase $3 of the gluino mass is irrelevant as well. Consequently, 
this associated production of chargino and neutralino and their subsequent decays involve 
only two CP-violating phases {$^,$1}. With these good approximations and after an 
appropriate Fierz transformation to the ul and dL-exchange amplitudes the transition 
matrix element can be written in the form 

T {du ^ xrx?) = -Q% imi.PMd)] [uix-)YPMx'j)] ■ (31) 



Here Qa/3 are the so-called generalized bilinear charges classified according to the 
chiralities a, f3 = L,R of the associated quark current and the chargino/neutralino current. 
The explicit forms of these bilinear charges are 

U - yy -Of^ 



ij _ 


Dw 




LL — 




ij _ 


Dw 




LR — 




ij _ 
RL — 


0, 


Qrr 



~ yy Df^ 
QlR = ^Rij + 

Q%L = 0, Q%R = 0. (32) 
with the S-, t-, and -u-channel propagators: 

Dw = — , D^ = -^, Dt^ = ^^, (33) 

s — + tmwi- w t — m^^ u — 
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where s = {pd + Paf, t = {pd - p^-f 



guj and gmj are given by 



and u = {pd —py,oY, and the couphngs WLij, Wffij, 



1 



Wmj = URiiNj2 - -^URi2Nj4, 



1 

V2' 



9Lij 



9Rij 



SwCw 

Uru 

SwCw 



sw{Qd + \)N*,-]^cwN*^ 



(34) 



Here Unj^Rij (Nij) are the chargino (neutrahno) mixing matrix elements and Qd = —1/3 
and Qu = +2/3 the electric charges for the down- and up-type quarks, respectively. Note 
that the coefficients guj and gjnj for the t- and w-channel diagrams are governed only by 
gaugino components of the chargino and neutralino while Wuj and are determined 

by both the gaugino and higgsino components. 



Defining the x7 flight direction with respect to the down quark momentum direction by 
B, the explicit form of the production helicity amplitudes can be determined from eq. 
In the limit of neglecting the initial d and u quark masses, the d and u helicities are opposite 
to each other in all the exchange amplitudes, but the x7 and helicities are less correlated 
due to the non-zero masses of the particles; the amplitudes with equal chargino/neutralino 
helicities oc m -- -o / ^/s must vanish only for asymptotic energies. Denoting the down quark 

A i ' Xj 

helicity by the ffist index, the x~ and x? helicities by the remaining two indices, the 



production helic: 
2-component sp 



+;++ 
+;+- 
+;-+ 
+; — 
-;++ 
-;+- 



ty amplitudes T{a; Aj, A^) = 27ra {a; XiXj) can be derived by the so-called 
nor technique 0], which yields 



Q'L^i-vl + Q'Ui-v 



sin 



Q'Mi^ + + V-) + Q'Uii - ^+)(i - V 



+ 



RR 



Qrr 
Qlr 



(1 - r^+)(l - r/_) + Q'^^^il + r^+)(l + r^_) 

sin 9 , 
sin 9 , 



^1 + cos 9) , 
;i -cos9). 



i-t- + Q'L\^i-n- 



i-vi + QlL\^i-v'- 



+ 



'1 — COS 9) , 
^1 + cos 9) , 



14 



-) = + 



sin 9 , 



(35) 



wliere r]± = A^/^(l, /if , /i^) ± (/if — /i^) witii /if (/i^) = m^„(m|o/s) and X{x,y,z) 

+ — 2xy — 2yz — 2zx. If tlie arguments are not specified, then the notation A stands 
for A(l,/if,/i^) in the following. 



All physical observables constructed through the production process are determined by 
the production helicity amplitudes and expressed in a very simple form by sixteen quartic 
charges [^] containing the full information on the dynamical properties of the production 
process. These quartic charges are expressed in terms of the bilinear charges given 
by eq. (^) and classified according to their transformation properties under parity (P) as 
follows: 



(a) Eight P— even terms: 



Qi 



1 r 

4 
1 

2 
1 

4 

- ^m 
2 

1 



^RR^RL "T ^LL^LR 



IQIlI' + \Q%r\' - \Q%l\' 



2 

1 

2 

3fJe 

3fJe 



Q'm 



^RR^RL "T ^LL^LR 
^RR^LR "T ^LL^RL 
^RR^LR ^LL^RL 



^rr'^ll 

^RL^LR 



\Q' 



LR\ 



(36) 



(b) Eight P— odd terms: 



Q2' 



3?e 



4 
1 

2 
1 

4 

- ^m 



\Qrr^ + IQiiP 



^RR^RL ^LL^LR 



^RR^RL ^LL^LR 



\Q1l\'] 



\Q1l\'] 
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53m 



^RR^LR ~ ^LL^RL 



53m 



^RR^LL 
^RL^LR 



(37) 



We note that these 16 quartic charges comprise the most complete set for any fermion-pair 
production process in du colhsions when the quark masses are neglected. On the other 
hand, the quartic charges defined by an imaginary part of the bilinear-charge correlations 
might be non-vanishing only when there exist complex CP-violating couplings or/and 
CP-preserving phases like rescattering phases or finite widths of the intermediate particles. 
Therefore, such nonvanishing values of those quartic charges may signal CP violation in a 
given process. 



4.3 Parton— level production cross section 

The unpolarized differential production cross section of the parton-level process du — > xiX^j 
is obtained straightforwardly by taking the average/sum over the initial/final helicities: 

'^^ {du - x-X?) = ^AV2 ^ |(^; A.A,)p . (38) 



dcosG' 32s 

Carrying out the sum, one finds the following expression for the parton level differential 
cross section in terms of the quartic charges: 

-^(rf^-xrx?) = ^AV2|[4-(r^+-r^_)' + (.^+ + r/_)2cos'e]g^ 
dcosH ■' 8s I 

+V(l-r/2)(l-r^2)g^. + 4^^^ + g^.j _ (39) 

Thus, the three P-even quartic charges Qi, Q2 and determine the 9 dependence of 
the cross section completely. 

The final cross section of the xix^j production in pp collisions is obtained by convoluting 
the parton level cross section with a parton distribution and it is given by, 

K dx r 1 

a{pp xlxl + ^) = — dr fd/p{x)fu/p{T/x) + fu/p{x)fd/pir/x) a{s = rs) ,(40) 



3 Jo X 

where / are the respective parton fiuxes in p and p and Tmin = {fri^- + m^of' / s. We have 
used the CTEQ4m [p3] parametrisation to obtain parton distribution setting the relevant 
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QCD scale to the cm. energy s of the parton level process. We take into account the dom- 
inant QCD radiative corrections to the production cross section by taking the parameter 
K = 1.3 in eq. (^). The production cross section for the associated positively-charged 
chargino and neutralino pair is the same as its charge-conjugate one. 

Figure || shows the distribution of the production cross section a{pp — > X1X2 + ^) with 
the scattering angle B for various sets of {^fi, $1} in the scenario SI and S2 for a given 
cm. energy of 1.8 TeV. We find that the production cross section is very sensitive to 
but (almost) insensitive to $1. In order to look into this feature more clearly, we present 
in Fig. the integrated production cross section (T{pp — > X1X2 + ^) ^^e {$^, $1} plane 
in (a) SI and (b) S2, for which the constraints from the electron EDM are embedded on 
the plane. The results exhibit a few interesting aspects: 

• The differential cross section is (almost) forward-backward symmetric in both scenar- 
ios. This can be understood by noting the fact that the higgsino parts of Wl^r are 
suppressed in both scenarios due to large values of and that the relevant quartic 
charges (5i,2 are forward-backward symmetric but Q3 is proportional to cos to a 
good approximation due to the assumption rriu = m^. 

• The large and small squark masses in the scenario S2 reduce the production 
cross section due to a severe destructive interference between the Pi^-exchange and 
squark-exchange diagrams. On the other hand, in the scenario SI, the t- and u- 
channel contributions with the large squark masses can be ignored so that the lack of 
destructive interference yields a much larger value of the production cross section in 
the scenario. 

• Except for the region of <l>^ = 0, 27r in SI, the integrated production cross section is 
almost independent of the phase $1 and it decreases as approaches vr in both sce- 
narios. This property is mainly set by the dependence of the chargino and neutralino 
masses on the CP-violating phases. 

Consequently, the production cross section itself can be a very sensitive probe to the phase 
but not to the phase $1. 

4.4 Chargino and neutralino polarization vectors 

The spin-1/2 chargino and neutralino through the parton-level process du X7X^ are in 
general polarized and their polarizations are reflected in the distributions of their decay 
products; x~ Xi^~^e and Xj ~^ Xi^'^^'~ ■ I^i order to have a rough estimate of the ef- 
fects of the CP-violating phases on the spin correlations, it is meaningful to investigate the 
chargino and neutralino polarizations in the parton-level production process du X7X^- 
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The polarization vector V^' = ("Pf^ ,V^' ,V^' ) of the produced chargino is de- 
fined in the rest frame in which the axis z\\L is in the fiight direction of x7 ^ ^\\T rotated 
counter-clockwise in the production plane, and y = z x x\\N of the decaying chargino x7 ■ 

Accordingly, V^' denotes the component parallel to the x7 flight direction in the cm. 

frame, V^" the transverse component in the production plane, and V^" the component 
normal to the production plane. The three components of the chargino polarization vector 
can be expressed by the production helicity amplitudes (|35| ) as 



= ^ E {K^; ++)!' + K^; +-)!' - K^; - K^; — /^^ 
'Pt' = E [(^; ++)(^; -+)* + — )(^; +-)*] }/-^^ 

'P¥ = 5: [{a; — ) {a; +-)* - {a; ++) {a; -+)*] j/AT , (41) 

with the normalization corresponding to the unpolarized distribution 

A^=^E [l(+;A.A,)r + |(-;A.A,)r] . (42) 

The polarization vector of the neutralino can be obtained similarly from the production 
helicity amplitudes (|3^) by exchanging the chargino and neutralino helicities in eq. (|4T|) 



'P? = \ E {K^; ++)!' + K^; - K^; +-)!' - K^; — )P} 

P^ = E [(^; — ) -+)* - ++) j/A/" • (43) 



The longitudinal and transverse components of the polarization vectors are P-odd and 
CP-even, but the normal component is P-even and CP-odd so that the normal polariza- 
tion component can be generated only by the complex production amplitudes. Certainly 
there exist non-trivial phases in CP non-invariant SUSY models. Also the non-zero width 
of the Z boson and loop corrections generate non-trivial phases; however, the Z boson 
width contribution to the normal polarization is negligible for high energies as mentioned 
before, and so are the radiative corrections. As a result the normal component is effectively 
generated by the genuine CP-odd phases of the couplings. 
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It is straightforward to find the normal polarization components of Xi and x^'- 

p¥ = ^ A^/V. sin e Q'i , 'Pn =^ A^/V, sin 9 Q'i . (44) 

As expected from the CP properties of the normal components, they are determined by the 
CP-odd quartic charge Q^^ , which requires some complex couplings. In order to estimate the 
size of the normal polarization quantitatively, we present in Fig. M the normal polarization 
Vj^ for four different combinations of {$^,$1} in the scenarios 51 and S2 taking the 
parton-level cm. energy of 300 GeV for the sake of illustration. [The chargino normal 
polarization is proportional to the neutralino normal polarization.] The expression ( PI ) 
and Fig. |] lead us to the following features: 

• The normal polarizations are suppressed near the thresholds or at high energies. 

• In the scenario SI with large squark masses, the quartic charge Q4 is forward- 
backward symmetric so that the normal polarizations also are forward-backward 
symmetric. However, this property is not maintained in the scenario S2. 

• In both scenarios, the size of the normal polarizations is too small to measure the CP- 
violating phases directly in the production of the associated chargino and neutralino 
pair. 

To conclude, it is likely that after applying the stringent EDM constraints to the CP- 
violating phases, we could not expect to have the normal polarizations of the chargino and 
neutralino large enough to be measured at the Tevatron. 



5 Polarized Chargino and Neutralino Decays 

The detection efficiency of the tri-lepton signatures at the Tevatron relies crucially on the 
branching fractions and distributions of the final three leptons, the latter of which depend 
on the polarizations of the decaying chargino and neutralino. To estimate the branching 
fractions, we need to calculate all the main partial decay widths that depend on the spar- 
ticle and Higgs-boson spectra in the MSSM. This section is devoted to a comprehensive 
discussion on the chargino and neutralino leptonic decays, including the polarizations of 
the decaying chargino and neutralino and the branching fractions. Firstly, we present the 
chargino and neutralino decay amplitudes in terms of the corresponding bilinear charges 
and the polarized decay distributions in terms of the quartic charges. Secondly, we calculate 
the decay density matrices by using the so-called Bouchiat-Michel formulas [^, that are 
to be used to form the complete production-decay spin/angular correlations. Finally, we 
estimate the branching fractions of the leptonic decays Xi ~^ x\^~^t and X2 ~^ Xi^'^^'~ 
in the scenarios 51 and S2. 
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5.1 Polarized decay distributions 



The diagrams contributing to the process x~ Xi^~^e. and Xj* 
/i are shown in Figs, ^(a) and (b). Here, the exchange of the neutral and charged Higgs 
bosons [replacing the W~ and Z bosons] are neglected since the Yukawa couplings to the 
light first and second generation leptons are very small. In this case, all the components of 
the decay matrix elements are, after a simple Fierz transformation, written as 



with i,e 



xie-i'+ 



(45) 



where a,(3 = L,R. Note that since the decaying neutralino is treated as an anti-particle 
in the associated production process du XiX% the v spinors for the decaying neutralino 
and the LSP appear in the expression for the neutralino leptonic decay. However, owing to 
the Majorana property of the neutralinos it does not matter whether the decay neutralino 
is treated as a particle or an anti-particle. The generalized bilinear charges C^^ for the 
chargino leptonic decay are given by 



LL 



C 



Di 



w 



RL 



V2i 



w 

^LR = CrR = 0) 



9Rily 



(46) 



The couplings Wui, Wru, Qlh, and qrh are given in eq. 
propagators are 



. The s'-, t'- and w'-channel 



w 



s' — rriyy + imw^w 



t' ~ m 



2 ' 



u' — m? 



(47) 



where the Mandelstam variables s\ t' and u' are defined in terms of the 4- momenta, go? Q 
and g, of x?, ^~ and Vi, respectively, as 



s' = {q + qY , t' = {qo + q) 



u 



{qo + q? 



(48) 



On the other hand, the generalized bilinear charges N-j^^ for the leptonic decay x 



X? d-''^^' are given by 



n" 



„2 2 



1 i' 
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Kr 



D'' 



' „2 ^jl ^t" "'Rjl ' 



where the s"-, t"- and ^''-channel propagators are 



D" — 



2 

mi 



^u" — 



U" 



2 



(49) 



(50) 



with the Mandelstam variables s" = {q' + q'y, t" = {q^ + g')^ and u" = [q^ + g')^ in terms of 
the 4-momenta, q'^, q' and q', of x?) and respectively. The couplings Zij, huj and 
hjiij are expressed in terms of neutralino diagonalization matrix elements N^j as 



■'ij 



huj 
hRij 



^ {Ni2Cw + Nasw) iN*cw + N*sw) , 



4S2 .2 

1 



N^lN*, 



and they satisfy the Hermiticity relations reflecting the CP relations 



2ij 



ji ' 



h 



Lij 



h* 

"'Lji • 



h 



Rij 



h* 

"'Rji ■ 



(51) 



(52) 



Note that Zji is governed by the higgsino components of x° while hiji and hRji arc deter- 
mined by the gaugino components of the neutralino. Therefore, as will be shown later, Z21 
are suppressed in the scenario S2 with a large while the t- and w-channel diagrams are 
suppressed in the scenario Si with large selectron masses. 



Applying the polarization projection operator of the chargino to the amplitude squared 
of the chargino leptonic decay Xi ~^ Xi^~^e yields the polarized decay distributions with 
the polarization vector nf of the chargino x~: 



\V\\n,) = -4(t'-m^_)(t'-m^o)(Ci-C^)-4(«'-m^-)(«'-mJo)(Ci + C^) 
—Sm^-m^os'Cl 

Xi Xi ^ 

-8(n, • q) [m--(m|o - u'){C^^ + Ct) + m^o(m|- - t')&^ 



+8(n, • q) [-m-- (m^o - t'){Ct - Ct) + m^o{ml- - u')Ct 



-16 m^o{qiniqq) CI , 



(53) 
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where {qiniqiq2) = ^^lupaQin'^q'^q'' with the convention 60123 = +1- Here, the quartic charges 
{Cji - Ci4} and {C^^ - C^^} for the chargino decays are defined by 



RR\ 



+ \cir.\' + \c],rr + \c 



RL\ 



CI 



CI 



-3fie 



'RR\ 



n'i- 1 2 

'^RU 



CI = -^Jm 
^ 2 



'^RR'^LR + 



LL'-^RL 



4 
1 

2 

'^3=7 



+ \c 



i |2 



- \C 



i |2 



i^i 1 2 



3f?e 



^RR^LR ^LL^RL 



\c 



RR\ 



\n'i- 1 2 
y^LR\ 



1 2 

y^RL\ 



\c 



LL\ 



(54) 



The polarized distribution with a polarization vector of the neutralino decay x° 
x5 can be derived in a straightforward way 

\V\\n^) = -4(i''-m|o)(t''-m|o)(7Vf + 7V|')-4K-^o)K-m|o)(7Vf-7V|) 



—8m^om^i)s"Nn 



-8{nj ■ q') 

+8(% • q') 

+ 16 m^o {qjnjq'q')Nl 



,.,,.40 - t"){Ni' + ATI') + ,n,o(m|o - u")N^^ 



■3 



(55) 



The quartic charges for the neutralino decay case can be obtained from the bilinear charges 
A^^^ in the same way as the chargino quartic charges are defined in terms of the bilinear 
charges C^^. Related to CP violation it is worthwhile to note that the quartic charges C4 
and Nl manifest CP violation in the theory. 

For the sake of subsequent discussion of the spin/angular correlations between the pro- 
duction and decay processes, wc construct the decay density matrix p\\i ~ VxDy. In 
general, the decay amplitude for a spin-1/2 particle and its complex conjugate can be 
expressed as 

V{X) = r u{q, A), V*{\') = uiq, A') f , (56) 

with the general spinor structure F and F = 7*'F''". Then we use the general formalism to 
calculate the decay density matrix involving a particle with 4-momentum q and mass m by 
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introducing three spacelike 4-vectors ra^ (a = 1,2,3) which together with q/m = form 
an orthonormal set: 
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ab 



9ab n^nl 



(57) 



with gf^'^, g"-^ = diag (1,-1,-1, —1) (a, 6 = 0, 1, 2, 3). A convenient choice for the exphcit 
form of n"^ is in a coordinate system where the direction of the three-momentum of the 
particle is g = (sin 9, 0, cos 9) lying on the x-z plane: 



n 



(O,cos^,0,-sine) 



n 



,2M 



(0,0,1,0) 



n 



m 



(5^ 



Then in the given reference frame, the three 4-vectors as defined in the above equation 
describe the transverse, normal and longitudinal polarization of the decaying particle, re- 
spectively. 



With the four-dimensional orthonormal basis of the 4-vectors {n ,n ,n ,n }, we can 
derive the so-called Bouchiat-Michel formulas p6|] for u and v spinors 

1 



u{q, X)u{q, A') = - [6xy + 75 ^r^^x] U + , 



v(q,\)v(q,X') = -[5,y+^, ^Xv] 



(59) 



with A, A' = ±. These formulas enable us to compute the squared, normalized decay density 
matrix pxy as follows: 



Paa' 



Vi\)V*{\') _ 1 



5xx' + 



'A'A 



(60) 



where r'^ (a = 1, 2, 3) are the Pauli matrices. The three functions X and (a = 1, 2, 3) for 
the chargino leptonic decay 

for the neutralino leptonic decay x° — Xi can be obtained easily as 



Xi^ and the three functions X and (a = 1,2,3) 



X 



-8(t' - ml-){t' - mJo)(q - CI) - 8{u' - mj_)(,.' - m^o)(Ci + C^) 
— \Qm--mc,os'Cl, 

-16« ■ q) [m^-(m|o - u'){Ct + Ct) + m^o(m^_ - i')Ct 

+ 16« • q) [-m-- (m|o - t')(C^f - ^sO + "^x;(<- - "O^s 
-32m^o(gi<gg)Q , 



(61) 
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and 



X = 



-8(i" - m|„)(t" - m|o)(iVf + Ni) - d>{u" - m|o)K - m|o)(7Vf - Ni) 
- 1 6 m^o m^o s"N2 , 



-16« ■ g') 



-m^o(m|o - u"){Ni' - Ni') + m^o(m|o - t")N'^ 
+m^o(m|o - t")(iVf + Ni!) + m^o(m|o - u")N'^ 



+32 m^o{qjn''q'q')Ni, 



(62) 



where n" and is the polarization vector of the decaying chargino and neutrahno, respec- 
tively. 



5.2 Branching fractions 

The main decay modes of the lightest chargino Xi and the next-lightest neutralino can 
be classified as follows: 



Xi 
XT 



W*xl H*xl^xlrv,, xiqq': 
B\ ut, qq' 

Z*x'i, H*xl 
£t,ui>*,qq*^xU^r,x'iqq, 



^i, Xiqq' , 
Xiqq, 



(63) 



with q and q' belonging to the same SU(2)i multiplet. Besides, if the mass m^± is smaller 
than the neutralino mass m^o, the lightest chargino xt can take part in the neutralino 
decay via the processes X2 ~^ Xi'^^*: Xi^^* ^^'^ vice versa. Concerning the main decay 
modes, there are several aspects to be noted: 



• For the first and second generation leptons, the Higgs-exchange diagrams are sup- 
pressed if tan P is not very large and there is no generational mixing in the slepton 
sector. 

• The experimental bounds on the Higgs particles are very stringent so that the two- 
body decays X2 ~^ Hxi, X2 ~^ H^x'f , and Xi ~^ H^Xa ^-i"© expected to be not 
available or at least strongly suppressed. If these decay modes open up, it could then 
spoil the tri-lepton signal. 

• The lightest chargino and the second-lightest neutralino are almost degenerate in the 
gaugino-dominated parameter space so that the charged decays such as X2 ~^ Xi '^^^e 
and Xi ~^ X2 ^~^e "W'ill be highly suppressed. 
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Keeping in mind the above subtle aspects, we calculate the branching fractions B{xi 
Xi(^~^i) and B{x2 ~^ x\^'^^'~) fully incorporating all the possible decay modes of the neu- 
tralino X2 but neglecting the Higgs-exchange contributions in both scenarios. 

We present in Fig. |10| the branching fractions B{xi Xi^~^i) and B{x2 ~^ Xi^'^^'~) 
for = e or /i in iSl (two upper figures) and in 52 (two lower figures). In the scenario 
SI the branching fraction B{xi — > Xi^~^^) is almost constant over the whole space of 
the phases and the branching fraction B{x2 ~^ x\^'^^'~) is very small and sensitive to the 
phases only around $1 around = 0, 27r. The insensitivity of both branching fractions to 
the phases in SI is due to the fact that the t- and w-channel contributions are suppressed 
due to large slepton masses and the couplings Wlii, Wmi, and Z21 for the s-channel con- 
tributions are not so much sensitive to and $1. On the contrary, the branching fractions 
are rather sensitive to $^ and $1 in the scenario S2. It is interesting that the branching 
fraction B{xi ~^ Xi ^~^e) can be minimal for certain non-trivial values of the CP-violating 
phases. We note that the branching fraction B{x2 Xi^'^^'~) is greatly enhanced in the 
scenario S2; this stems from the fact that the slepton-exchange contributions due to mainly 
the gaugino components of the neutralinos become dominant for the small slepton masses 
and the large value of On the contrary, the branching fraction B{xi — ^ x\^~^t) is not 
so different in size between the two scenarios, but the branching fraction becomes more 
sensitive to the CP-violating phases in the scenario 52. 

In summary, the branching fractions B{xi ~^ Xi^~^e) and B{x2 ~^ Xi^'~^^'~) are not 
so strongly dependent on the CP-violating phases {$^,$1}, but the branching fraction 
B{x2 ~^ Xi^'^^'~) is very sensitive to the slepton masses. 



6 Spin/ Angular Correlated Observables 
6.1 Correlations between production and decay 

In this section we provide a general formalism to describe the spin/angular correlations 
between the production process du Xix'j and the sequential leptonic decays of xi and 
Formally we can have the spin/angular correlated distribution by taking the sum over 
the helicity indices of the intermediate chargino and neutralino states and folding with the 
chargino-neutralino leptonic decay density matrix pw and pxy with the matrix squared 
for the production helicity amplitudes: 



E 



2 2 

TT a 



EEE(a;AA)(a;A'A')WpAA-' 



AA' AA' 



2 2 
TT a 



PzP + + + PyU + P.,V + PyV + P.P. Q 
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+ P,P,W + PyP.W + P^P.X + PyP.X + (P,P, - P,Pj,)3^ 

+(P,.P, + P.P.) 3^ + (P,.P, + P,P,) Z + (P.P.. - P,.P,) Z 



(64) 



where the functions Px,y,z and Px,y,z depending on the chargino and neutrahno decay dis- 
tributions, respectively, 



P. 




Py = 


Y2 

X ' 


P. 




p. 


~ X ' 


Py = 


X ' 


P. 


_ 



(65) 



serve as the polarimeters to extract the spin-spin correlations of the chargino Xa and 
the neutralino Xj* in the production process. The sixteen coefficients are combinations of 
the production lielicity amplitudes, corresponding to the unpolarized cross section, 2x3 
polarization components and 3x3 spin-spin correlations. 



I^i) Unpolarized part: 



-'unpol 



Ka;++)p+Ka;+-)p + |(a;-+)p + |(a; 
(ii) Polarization components: 



(66) 



V 
V 

u 

V 



E 



^ (T=± 

-E 



\ E 5Re{(a; ++)* + (a; — )(a; +-)*} , 

i 3?e|(a; ++)* + (a; — -+)*} , 



(67) 



and W, V defined as W, V after replacing 3fJe by 55m. We note in passing that the above com- 
binations are directly related with the polarization vector of the chargino X\ and neutralino 
X2 defined in Sect. ^|4| 



^;-+ r+ ^x: 



(iii) Spin— spin correlations: 

1 



4 



E 

cr=± 
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^ = ^E5Je{(^;-+)(^;++)*-(^;— )(^;+-r}, 
^ = ^E^(^;+-)(<^;++>*-(<^;—>(<^; 

i^^5Re{(a;-+)(a; +-)*}, (68) 



Z - 

and VV,A',3),Z defined as W , X ,y , Z after replacing 3?e by S^m, and these components 
along with U and V will contribute to CP-odd observables. 



Combining the production and decay distributions, we obtain the fully spin/angular 
correlated 11-fold differential cross section for the parton-level process du —>■ xiX^j ~^ 

d<7 = ^ B{x- - Xl r P,) B{x] - X? e-e^) E d$3. , (69) 

°^ corr 

where d$3£ denotes the final-state phase space volume element and it can be parameterized 
in terms of 11 independent kinematical variables as follows: 

d$3^ = d cos © dxi dx2 d cos 9i d^i d0i2 dxs dx4 d cos 6^ d^a d034 , 

where the angular variable 9i is the polar angle of the ^'' in the xi rest frame with respect 
to the original flight direction in the parton-level du center of mass frame, and 01 the 
corresponding azimuthal angle with respect to the production plane, and 0i2 is the relative 
azimuthal angle of Vf along the £^ direction with respect to the production plane. A similar 
configuration can be specified for the neutralino decay distribution by ^^3, 03 and 034. The 
dimensionless parameters xi, X2, 0:3, and denote the lepton energy fractions 

Xi = , X2 = , X3 = , Xi = . (70) 

m^- m^- m^o m^o 

Xi Xi Xj Xj 

The allowed space of the kinematical variables for the chargino leptonic decay is determined 
by the kinematic conditions obtained with the masses of the final-state leptons neglected: 

< e < tt; < ^1 < tt , < 01 < 27r , < 0i2 < 27r; 
0<xi,2<l-ra, (1 - xi){l - X2) > ra , Xi + X2> I -ra, (71) 

where rn = m|o/m?_, and similarly the allowed range of the kinematical observables for 
the neutralino leptonic decay can be obtained simply by replacing the labels; (1, 2) to (3, 4) 
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and rn r,i = m?o/m|o. 

Finally, the tri-lepton rates in the pp laboratory frame is obtained by folding the parton- 
level differential cross section (^) with the d-quark and w-quark parton distributions. This 
folding involves 2 additional kinematical variables r and x. Since the parton-level cm. 
frame is not fixed with respect to the pp cm. frame the Lorentz boost of the partonic 
system along the initial beam direction has to be properly taken into account. 



6.2 Total cross section of the correlated process 

The total cross section for the correlated process pp —>■ 31 + X can be obtained simply by 
computing 

cripp ^3i + X) = a{pp ^ ^7 - X? ^ z7,) B{xl - xl ^'^n • (72) 

We present in Fig. [Tl| the contour plots for the total cross section on the plane of two 
CP-violating phases {$^, $1} in the two scenarios (a) SI and (b) 52. In the scenario 51 the 
branching fraction B{xi ~^ Xi is almost constant as shown in Fig. 10 so that the total 
tri-lepton cross section is mainly determined by i3(x2 ~^ Xi ^'~^^'~) and the production cross 
section a{pp — > xFX^)- The total tri-lepton cross section is of the order of 1 fb (10 fb) for the 
parameter sets 51 (52), respectively. However in the future Tevatron RUN-II experiments 
with the upgraded luminosity of the order of 2 fb~^ a handful of tri-lepton events can 
be produced. Note that the total cross section is sensitive to the CP-violating phases in 
the scenario 52 and the cross section can be minimal for certain non-trivial values of the 
CP-violating phases. This property is mainly due to the fact that the branching fraction 
B{xi Xi^~^i) exhibits a similar pattern as shown in Fig. [TO. Therefore, depending on 



the size of the integrated luminosity the very existence of the minimum event rate and 
the simultaneous small mass splitting (See Fig. ^ for the non-trivial CP-violating phases 
reflect that the Tevatron bounds on m^± and m^o might be much smaller than those p7 
ruled out in the context of SUGRA and GUT inspired SUSY models. 



6.3 Dilepton invariant mass distributions 

The final-state leptons of the neutralino decay X2 ~^ Xi^'^^'~ provides us with a very 
easily measurable kinematical observable; the dilepton invariant mass, mu. This Lorentz- 
invariant quantity can be precisely reconstructed by measuring the two lepton momenta, 
and it is nothing but the square root of the Mandelstam variable, a/s", 

rriu = \fs^' = m^o^xs + X4 — 1 + rji . (73) 

Furthermore, the invariant mass distribution is independent of the specific production pro- 
cess for the parent neutralino X2^ because the invariant mass does not involve any angular 
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variables describing the decays so that the polarization of the decaying neutralino does not 
affect the distribution. 



Figure |T2| shows the dilepton invariant mass distribution in the scenarios (a) SI and (b) 
iS2. In principle, three dilepton combinations can be constructed out of the three charged 
leptons. But, the correct dilepton combination will have a sharp end point of its invariant 
mass distribution with a distinguishable peak in the dilepton invariant mass distribution. 
This allows one to experimentally determine the mass difference = m^o — m^o with a 
good precision. Note that the position of the end points is strongly dependent on the CP- 
violating phases. It can provide us with the opportunity to probe the CP violating phases 
if all the other real parameters are known. As the gaugino masses vary significantly with 
the relevant CP violating phases, the distributions with respect to the kinematic variables 
like the energies and transverse momenta of the final-state leptons are strongly influenced 
by those CP-violating phases. 



6.4 Lepton angular distribution in the laboratory frame 

In a realistic experimental situation, it is necessary to isolate tri-lepton signal events by 
applying various selection cuts effectively to reduce the contaminations from all the back- 
ground processes as much as possible. For that purpose, it is very important to fully 
understand the event topology of tri-lepton signal which depends on the polarizations of 
the chargino and neutralino at the intermediate stage. The full incorporation of the spin- 



correlations (|6^) involves a lot of correlated terms, which many Monte Carlo simulations 
T3[] have simply neglected by including only the first un-correlated term in eq. (Q). In 
this section, we take an easily-measurable kinematical variable, the scattering angle 6i of 
the lepton from the chargino decay with respect to the proton beam direction and estimate 
the variation of the lepton angular distribution of the tri-lepton signal due to the spin- 
correlation effects. 



Figure 13 exhibits the lepton angular distribution of the tri-lepton signature for both 
the un-correlated case and the fully-correlated case for different phase combinations in both 
scenarios. Note that the correlated distribution can be different from the the un-correlated 
distribution around the forward and backward regions depending on the values of the CP- 
violating phases. In particular, the forward-backward asymmetry can be different in two 
cases. Therefore, even for this simple distribution it is important to take into account the 
spin correlations between production and decay fully to obtain the magnitude properly. 
Without any experimental cuts, such observables as the total tri-lepton event rate can 
be independent of whether the spin/angular correlations are taken into account or not. 
However, it might be inevitable to apply some efficient experimental cuts to suppress serious 
backgrounds. In that case the spin-correlations should be included. 
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6.5 CP— odd triple momentum products 



So far we have concentrated mainly on the CP-even production-decay correlated observ- 
ables which depend on the CP phases only indirectly. For direct measurements of the 
CP-violating phases one has to use CP-odd or T-odd observables. Some of such CP-odd 
(or T-odd) observables can be constructed by taking a triple product of any combination 
of the initial proton (or anti-proton) momentum and the three final lepton momenta. One 
typical example is the following triple momentum product (TMP): 

Ot=P,^-{p,^ xp;j, (74) 

where ii = i~ of the chargino decay xF ~^ x\^~^i^ and £3 = £4 = of the neu- 
tralino decay X2 ~^ Xi^'~^'^- The observable (^) enables us to probe the CP-violating 
phases directly when neglecting the tiny particle decay widths. Similarly, the initial proton 
momentum and two final-state leptons allows us to construct additional T-odd observables: 

Ot = Pp ■ {Pe X Pi') (75) 

where {i, i'} is any combination of two momenta among the three final lepton momenta. 
In total, we have four independent TMPs; Ot, Ct'^', Ct'^' and Ot'^\ 



In general, any of T-odd TMP can be given by a linear combination of the quar- 
tic charges {Q4,Cl,Nl}. We note that the same topological pattern of the contributing 
diagrams between the associated production and the chargino decay make and C| cor- 
related in size. Even before making any numerical estimate of the T-odd triple products, 
we can argue that their size is very small in the scenario SI with heavy sfermion masses. 
Firstly, the chargino and neutralino normal polarizations are very small in the scenario 
SI as shown in Sect. ^I4| , implying small and C4. Secondly, with the negligible t- 
and M-channel slepton contributions, the remaining quartic charge N4 of the neutralino X2 
simplifies to the expression 



\D 



^ IsL--] -sm 



{Zli) , (76) 



which contains a very small numerical factor (s^ — 1/4) ~ —0.02 PB[. Therefore, the 



quartic charge A^4 is also extremely suppressed in the scenario SI. This simultaneous sup- 
pression of three quartic charges render all the TMPs strongly suppressed in the scenario SI. 

The three TMPS {Ot, Ot^^, Ot^^} involve both the chargino and neutralino leptonic 
decays, but the TMP O^''^^^ involves only the neutralino leptonic decay. Therefore, one 
can have a large statistical gain by exploiting the observable Ot^"^ in measuring the CP- 
violating phases. In this light, we consider only the observable Oj?^" to probe the CP-odd 
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phases which are not excluded by the EDM constraints. Certainly, one needs to estimate 
all the possible systematic uncertainties in determining the reconstruction efficiency of the 
X1X2 mode. Nevertheless, let us take into account only the statistical errors including the 
hadronic decay modes of the chargino the present work in which case the excluded 

region of {$^,$1} at the N-a level for a given integrated luminosity J Cdt satisfies the 
inequality: 



where atot = a{pp ^ xFx") ^(x" ^ and {X) = / X ^d^/a^i over the total 

phase space volume $. The numerical factor 2 in the denominator is due to two possible 
combinations of the two final-state leptons; (e~,e+) and 



Figure exhibits the region of the CP-violating phases and $1 that is excluded 
by both the electron EDM constraints at 95% confidence level (shaded region) and by the 
T-odd observable O^'" at the 2-cr level with an integrated luminosity of 20 fb~^ (filled 
circles) and 30 fb~^ (open circles) in the scenario S2. Certainly, after incorporating all the 
systematic errors, the covered region should be reduced. Nevertheless, it might be useful to 
measure the T-odd observable at the upgraded Tevatron with an integrated luminosity C 
of the order of 30 fb~^ because the electron EDM and the T-odd TMP are complementary 
in constraining the CP-violating phases. 



7 Conclusions 

In this paper, we have investigated in detail the impact of the phases and $1 on the 
SUSY tri-lepton signals at the Tevatron in the framework of MSSM with general CP phases 
but without generational mixing. The stringent constraints by the electron and neutron 
EDM on the CP phases have been also included in the discussion of the effects of the CP 
phases. For the sake of illustration, we have considered two exemplary scenarios for the 
relevant SUSY parameters; 51 with very heavy first- and second-generation sfermions and 
iS2 with relatively light sfermions but a large 

We have found that in both scenarios the CP-violating phases can have a significant 
impact on the production cross section and the partial leptonic branching fractions of the 
chargino xf and neutralino X2- As a result, there may lead to a minimum rate of the 
tri-lepton signal for non-trivial CP phases. This implies that one should be careful when 
interpreting the chargino and neutralino mass limits derived under the assumption of van- 
ishing phases, since the worst case is not (always) covered by just flipping the sign of fi; 
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rather it can occur from some non-trivial phases in between. 

The production-decay spin correlations lead to several CP-even observables. We have 
studied the useful kinematical observables such as the dilepton invariant mass distribution, 
the angular distributions of the final-state leptons, and the T-odd (CP-odd) triple mo- 
mentum products of the initial proton momentum and two final lepton momenta. 

We have found that the the end point of the dilepton invariant mass distribution is 
very sensitive to the relevant CP-violating phases because of the strong dependence of 
the neutralino masses on the phases. Therefore, these distributions can be very useful in 
determining the phases once the other real SUSY parameters are known. The angular dis- 
tributions of the final-state leptons taking into account the full spin/angular correlations 
can differ from the non-correlated ones by a few percents. Therefore, it will be sometimes 
necessary to consider the fully-correlated distributions in order to interpret experimental 
data properly. 

It turned out to be difficult to investigate the CP-violating phases directly through the 
T-odd triple momentum products at the Tevatron with its upgraded luminosity of about 
2 fb~^. But, we have found that a substantial region of the CP-violating phases may be 
explored through the triple momentum products with the luminosity of about 30 fb~^ as 
proposed for TeV33. 
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(a) (b) 



Figure 1: The Feynman diagrams contributing to the fermion EDMs and CEDMs in the 
MSSM; (a) one-loop chargino, neutralino and gluino diagrams and (h) two-loop top-Higgs 
or top quark-top squark-gluino diagrams. 




Figure 2: (a) the allowed range of the CP-violating phase versus the higgsino mass 
parameter and (b) the allowed region ci/{$^, $1} against the electron (filled circles) and 
neutron (hollow circles) EDM constraints in the scenario S2. The trilin ear parameter \ Ae\ 
is taken to be 1 TeV and its phase is scanned over the full allowed range. 
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(a) x"i Mass [GeV] (b) % ^ Mass [GeV] (c) % ^ Mass [GeV] 



^ _|_ 

40 ► >*: 



+ 70 
X 85 
• 100 



+ ''x * 

+ y. % 

+ X • 

+ \ • 

+++ ^- • 



V 



X 45 
• 50 



S \ 



Xx > 
X % 

^ i 

+++ 

V + + 
\ + 

X + 

k << +, 



+ 70 
X 85 
• 100 



95 
100 
i 104 



+ 49 , "U "x 4, 
>-50 \ 
'51 \\\ \ 

V >>< •, 



+ 95 
X 100 
• 104 



• >>: + 



+ \ 
+ 

+ xj 

+ 



• X +^ 

* X + 



60 



2 4) 4 



60 



2 ^> ^ 
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6A =400MeV 

m 




Figure 4: The constraints on {^ix, $1} by the measurements of the mass difference = 
m^o — m^o for four different mean values of the mass difference with the assumed 
uncertainty of 4OO MeV in (a) Si and (b) S2. 




Figure 5: Three mechanisms contributing to the parton-level production process du 
XiX^j', the s-channel W~ exchange, the t-channel ul exchange and the u-channel di ex- 
change. 
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(a) Scenario 1 



(b) Scenario 2 



2000 




1000 



© 

o 
u 

5 



■o 




Figure 6: The differential production cross section do'prod(pp ^ Xi X2 + -^)/dcos6' with 
respect to the cosine of the scattering angle, cos 9, in (a) SI and (b) S2 for three dijferent 
sets o/{$^,$i}. 
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(a) Scenario 1 



(b) Scenario 2 
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Figure 7: The production cross section a{pp 
the scenario SI and (b) the scenario S2. 



Xi X2 + ^) on the {^n, $1} plane in 
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(a) Scenario 1 



(b) Scenario 2 




Figure 8: The normal polarization component of the neutralino X2 ^^^^ respect to the cosine 

of the scattering angle, cos 0, in the process du xrX2 for four different combinations of 
{$^,$1} in (a) S\ and (h) S2. The parton-level cm. energy is assumed to be 300 GeV 
for the sake of illustration. 
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Figure 9: The three mechanisms contributing to (a) the chargino decay xi Xi^~^£ 
(b) the neutralino decay x° x^^'+f'". In both decays the exchanges of the charged and 
neutral Higgs bosons are neglected because they involve tiny lepton Yukawa couplings. 
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Figure 10: The branching fractions B{xi X? ^ ^e) B{x2 ~^ Xi ) for 
or n on the {$^, $1} plane in the scenarios SI (upper part) and S2 (lower part). 
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(a) Scenario 1 



(b) Scenario 2 
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Figure 11: The total cross section a{pp 3^ + X) of the tri-lepton signal on the {^^, $i 
plane in the scenarios (a) SI and (b) S2. 
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Figure 12: The differential correlated cross section datot/dmu with respect to the invariant 
mass rriii of two leptons from the neutralino decay X2 — > x1 for three different sets of 

{$^,$1} in the scenarios (a) SI and (b) S2. 
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Figure 13: The angular distribution of the charged lepton from the lightest chargino in 
the tri-lepton signal for the phases {^fj,, {0, 0}, {7r/2, 0}, {0, 7r/2} and {tt, 7r/2} in the 
scenarios Si and S2. 
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Figure 14: T/ie region of the CP-violating phases {$^,$1} excluded by the electron EDM 

measurements (black shadowed region) and probed by the T-odd TMP 0^^" measurements 
with the integrated luminosity of 20 fb~^ (filled circles) and 30 fb~^ (open circles) at the 
2-a level in the scenario S2. 
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